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Abstract. The small-scale power spectrum of the cosmological matter distribution, 
together with other cosmological data, provides a sensitive measure of the hot dark 
matter fraction, leading to restrictive neutrino mass limits. We extend this argument 
to generic cases of low-mass thermal relics. We vary the cosmic epoch of thermal 
decoupling, the radiation content of the universe, and the new particle's spin degrees 
of freedom. Our treatment covers various scenarios of active plus sterile neutrinos or 
axion-like particles. For three degenerate massive neutrinos, we reproduce the well- 
known limit of m v < 0.34 eV. In a 3+1 scenario of 3 massless and 1 fully thermalized 
sterile neutrino we find m„ < 1.0 eV. Thermally produced QCD axions must obey 
m a < 3.0 eV, superseding limits from a direct telescope search, but leaving room for 
solar eV-mass axions to be discovered by the CAST experiment. 
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1. Introduction 

The observed universe is surprisingly well described by a simple Friedmann-Robertson- 
Walker model that is characterized by a handful of well-measured global parameters. 
Together with the gravitational instability theory for the formation of structure, a 
"concordance model" has emerged that does not seem to conflict with any firmly 
established observations. Of course, the physical nature of most of the gravitating 
matter and energy remains elusive, except that most of the dark matter must be of the 
cold variety to avoid excessive free streaming in the early universe that would suppress 
small-scale power of the observed matter distribution. In the framework of this standard 
paradigm, the observed power spectrum of the cosmic matter distribution provides 
restrictive limits on the hot dark matter fraction of the cosmological matter cocktail [1]. 

Taking advantage of the latest cosmological data, notably the 2dF and SDSS galaxy 
redshift surveys and the WMAP measurement of the cosmic microwave background 
(CMB) temperature fluctuations, several authors have used this argument to derive 
restrictive neutrino mass limits of J2 m u < 0.7-1.0 eV (95% CL) [2-8]. The exact value 
of the nominal limit depends on the included data sets and on the assumed priors 
for some of the cosmological parameters. However, the small spread between different 
nominal limits suggests that the universe is a robust laboratory for constraining the 
overall neutrino mass scale, and perhaps eventually for positively determining its value 
[9-12]. Of course, these limits or a future positive detection depend on the validity of 
the standard cosmological paradigm. For example, the spectrum of primordial density 
fluctuations may not follow a simple power law, modifying all of these results. Still, 
very recently it was shown that the neutrino mass limit is surprisingly robust against 
a contribution of nonadiabatic, incoherent fluctuations such as would be predicted by 
topological defects [13]. 

We are here concerned with a different extension of the standard picture in that 
we study the effect of a non-standard number density of neutrinos, perhaps caused by 
the existence of additional neutrinos such as sterile states, or new low-mass particles 
such as axions. One of us has previously shown that increasing the assumed number 
density of neutrinos actually weakens the mass limit, contrary to naive intuition, because 
the increased radiation density counteracts the effect of the increased hot dark matter 
fraction [4]. Therefore, the usual neutrino mass limits do not trivially translate into 
more general cases where the hot dark matter fraction and the radiation content are 
varied independently. 

The free-streaming effect of low-mass particles on the power spectrum of the matter 
distribution and of the CMB temperature fluctuations depends on the particle mass mx, 
their number density and velocity distribution, and on the radiation content in other 
degrees of freedom that do not carry mass. We will consider two generic classes of cases. 
One covers neutrinos, i.e. fermionic degrees of freedom with a momentum distribution 
given by that of ordinary neutrinos. We vary the mass m u , the effective number of 
flavors carrying that mass, and the effective number of massless flavors. Our second 
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generic case is that of particles that were once in thermal equilibrium (thermal relics) so 
that their contribution to the hot dark matter fraction and their velocity distribution is 
perfectly defined by the cosmic temperature Td when they thermally decoupled from the 
background medium, by their mass, and by their number of internal degrees of freedom. 

In Sec. 2 we describe our general methodology for deriving the allowed range of non- 
standard cosmological parameters. In Sec. 3 we apply our method to various scenarios 
involving active and sterile neutrinos. In Sec. 4 we consider general thermal relics, and 
finally we conclude in Sec. 5. 

2. Likelihood Analysis 

2.1. Theoretical Predictions 

In order to derive limits on the parameters characterizing the new particles we compare 
theoretical power spectra for the matter distribution and the CMB temperature 
fluctuations with observational data in analogy to the previous works by one of us [3, 4]. 
The predicted power spectra are calculated with the publicly available CMBFAST 
package [14]. This package allows one to include directly the effective number of massless 
neutrino degrees of freedom and of neutrinos carrying a mass as well as the cosmic hot 
dark matter fraction, assumed to be equally distributed among the massive flavors. Put 
another way, the standard version of the code assumes a velocity distribution of the 
neutrinos that ordinary active neutrinos with a mass would have. For our more general 
case of arbitrary thermal relics we modified the code accordingly. 

As a set of cosmological parameters apart from the new particle characteristics 
we choose the matter density Q m , the baryon density Qf,, the Hubble parameter H , 
the scalar spectral index of the primordial fluctuation spectrum n s , the optical depth 
to reionization r, the normalization of the CMB power spectrum Q, and the bias 
parameter b. We restrict our analysis to geometrically flat models Q = Q m + Q\ = 1. 
Some of the other parameters are not kept entirely free, but are varied over a prior 
range that is determined from cosmological observations other than CMB and LSS. In 
flat models the matter density is restricted by observations of Type la supernovae to 
be fl m = 0.28 ± 0.14 [15]. Furthermore, the HST Hubble key project has measured 
Hq = 72 ± 8 km s _1 Mpc~ x [16]. The cosmological parameters and their assumed 
priors are listed in Table 1. The actual marginalization over parameters was performed 
using a simulated annealing procedure [17]. The cosmological parameters in our model 
corresponds to the simplest ACDM model which fits present data. While our results 
would be changed slightly by including additional parameters, there would be no 
qualitative changes. The neutrino mass is mainly degenerate with the matter density, 
the bias parameter, and to a lesser extent with the Hubble parameter, but not for 
instance with the curvature [3]. However, our prior on Q m from SN la data is not very 
restrictive, and the results would not change significantly by removing it. Removing the 
prior on H also will not lead to any substantial change in our results. 
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Table 1. Priors on cosmological parameters used in the likelihood analysis. 



Parameter 


Prior 


Distribution 


ft = fl m + Oa 


1 


Fixed 


Cirri 


0.28 ±0.14 


Gaussian 


h 


0.72 ±0.08 


Gaussian 


n b h 2 


0.014-0.040 


Top hat 


n s 


0.6-1.4 


Top hat 


T 


0-1 


Top hat 


Q 




Free 


b 




Free 



2.2. Cosmological Data 

2.2.1. Large Scale Structure (LSS). At present there are two large galaxy surveys of 
comparable size, the Sloan Digital Sky Survey (SDSS) [18, 19] and the 2dFGRS (2 degree 
Field Galaxy Redshift Survey) [20]. Once the SDSS is completed in 2005 it will be 
significantly larger and more accurate than the 2dFGRS. At present the two surveys 
are, however, comparable in precision and in the present paper we use data from the 
2dFGRS alone. 

Tegmark, Hamilton and Xu [21] have calculated a power spectrum, P(k), from this 
data, which we use in the present work. The 2dFGRS data extends to very small 
scales where there are large effects of non-linearity. Since we calculate only linear 
power spectra, we follow standard procedures and use only data on scales larger than 
k = 0.2/iMpc -1 , where effects of non-linearity should be minimal (see for instance 
Ref. [18] for a discussion). With this cut the number of data points for the power 
spectrum reduces to 18. 

2.2.2. Cosmic Microwave Background. The temperature fluctuations are conveniently 
described in terms of the spherical harmonics power spectrum Ct,i = {\aim\ 2 }, where 
^r(9,(j)) = J2im a imYi m (9 , <fi) . Since Thomson scattering polarizes light, there are also 
power spectra coming from the polarization. The polarization can be divided into a 
curl-free (E) and a curl (B) component, yielding four independent power spectra: Ctj, 
Ce,u Cb,u an d the T-E cross-correlation C T e,i- 

The WMAP experiment has reported data only on C T ,i and C T e,i as described in 
Refs. [6, 22-26]. We have performed our likelihood analysis using the prescription given 
by the WMAP collaboration [6,22-26] which includes the correlation between different 
C/'s. Foreground contamination has already been subtracted from their published data. 

In addition we use other CMB data from the compilation by Wang et al. [27] which 
includes data at higher I. Altogether this data set has 28 points. 
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3. Neutrinos 

3.1. Mass Equally Distributed Among N v Species 

Our first class of cases consists of different scenarios with N v flavors of neutrinos with 
equal momentum distributions at those early epochs where they are relativistic. Several 
of these species are assumed to carry equal masses m v . 

The first and simplest example is that of N v flavors carrying equal masses m v . In 
this case the hot dark matter fraction is given by the standard expression 

where h is the Hubble constant in units of 100 km s _1 Mpc -1 . Of course, the limiting 
case N v = 3 corresponds to ordinary neutrinos with a degenerate pattern of masses. 

A value for N v much larger than 3 contradicts the observed primordial helium 
abundance in the framework of the standard big-bang nucleosynthesis (BBN) theory. 
However, if the primordial v e chemical potential takes on a suitable value, its effect 
on the n/p freeze-out abundance can always compensate for the expansion-rate effect 
caused by the additional flavors. Such a fine-tuned v e degeneracy parameter may not 
seem particularly plausible, but it serves as one possibility to escape the BBN limits. In 
any event, we think it is useful to study what can be learned about neutrino parameters 
from large-scale structure data alone, without combining this information with data 
pertaining to an entirely different cosmological epoch. 

Since we know that the number of ordinary neutrino flavors is exactly 3, the 
additional degrees of freedom would have to be sterile states, perhaps Dirac partners of 
the ordinary neutrinos that were thermally equilibrated by some unknown mechanism. 

Without sterile states, in principle the effective N v can be enhanced by suitable 
neutrino chemical potentials for the ordinary flavors. However, the observed neutrino 
mixing parameters imply that all three standard flavors reach full thermal equilibrium 
(i.e. kinetic and chemical equilibrium) before the BBN epoch so that the restrictive BBN 
limit on the v e degeneracy parameter applies to all flavors [28]. Therefore, interpreting 
N v > 3 in terms of an increased density of the ordinary flavors requires some unspecified 
way to escape the BBN limits. We also note that our treatment assumes a neutrino 
velocity distribution corresponding to a non-degenerate thermal distribution. Therefore, 
a case with large neutrino chemical potentials is not strictly covered here. 

After marginalizing over the cosmological parameters shown in Table 1, the 68% 
and 95% CL allowed region of neutrino parameters is shown in Fig. 1. In the left panel 
this allowed region is shown in the N v -Q v h 2 -p\ane, in the right panel the equivalent 
region in the A^-m^-plane. 

From the right panel of Fig. 1 it is apparent that the 95% CL upper limit on m u 
is nearly independent of N u , contrary to naive intuition. One of us first observed this 
effect in Ref. [4]. The limits stated there agree with our present ones that are given 
explicitly in Table 2 for several values of N v . 
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Figure 1. Likelihood contours (68% and 95%) for the case of N v neutrinos with equal 
masses m v . Left panel: ./V^-f^/^-plane. Right panel: Equivalent -ZV„-m„-plane. 



Table 2. Neutrino mass limits (95% CL) for N u flavors with equal mass m, 
Flavors Mass limits [eV] 



N v 


N v m v 




3 


1.01 


0.34 


4 


1.38 


0.35 


5 


2.12 


0.42 


G 


2.69 


0.45 



If we specifically assume that the neutrinos are Dirac particles with a degenerate 
mass spectrum, and that the right-handed partners were thermally equilibrated by some 
unknown mechanism so that they have equal number densities with the left-handed 
states, the cosmological mass limit on m v degrades from 0.34 eV to 0.45 eV. 



3.2. Mass in One Fully Thermalized Species 

As a next case we consider one fully excited neutrino flavor that carries a mass m v so 
that the cosmic hot dark matter fraction is given by 

In addition there are N v — 1 massless species that contribute to the radiation density. The 
small mass differences between the ordinary neutrinos imply that this scenario cannot 
be realized among the ordinary neutrinos so that we need additional sterile states. For 
example, the 3+1 model for explaining the current neutrino oscillation data including 
LSND would be a case in point. As in the previous case, one must assume that the 
BBN limits are circumvented by some mechanism, for example a suitable v e chemical 
potential. Since we have 3 flavors of ordinary neutrinos that would contribute to the 
tally of massless states, and because one fully thermalized species is assumed to carry 
a mass, the present scenario requires N v > 4. The allowed region of iVVm„-values is 
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Figure 2. Likelihood contours (68% and 95%) for the case of N v flavors with exactly 
one of them carrying a mass m„. Left panel: N l/ -Q,^h 2 -plane. Right panel: Equivalent 
^-TO^-plane. 

Table 3. Neutrino mass limits (95% CL) for N v flavors with exactly one of them 
carrying mass m„. 



Flavors 


Mass limit [eV] 


3 


0.73 


4 


1.05 


5 


2.47 


6 


4.13 



shown in Fig. 2 where we have extended the plot down to N v = 3. In Table 3 we give 
explicit numerical 95% CL mass limits for several values of N v . 



3.3. Three "Massless" Standard Neutrinos 

Our last neutrino example is one where we consider the three standard neutrinos to 
have very small hierarchical masses, i.e. to be effectively massless for our purpose. In 
addition we assume there to be N v — 3 other flavors, e.g. sterile states, that carry equal 
masses m u so that the hot dark matter fraction is 

n„h 2 = (N„-3) m \ . (3) 

V ; 91.5 eV V ; 

An example would be a scenario with an additional sterile neutrino that did not reach 
full thermal equilibrium so that it is present with less than one full effective degree of 
freedom. 

The allowed region of N v -m v - values is shown in Fig. 3. In Table 4 we give explicit 
numerical 95% CL mass limits for several values of N u . Note that the case N v = 4 is 
identical to the N v = 4 case in the preceding section. The slight difference in quoted 
mass limit (1.02 eV instead of 1.05 eV) is due to the fact that the limit is interpolated 
from the calculated likelihood grid. 
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Figure 3. Likelihood contours (68% and 95%) for the case of N v flavors, N v — 3 each 
carrying mass ro„, and an additional 3 massless flavors. 

Table 4. Neutrino mass limits for the case of N v species, with 3 massless flavors and 
N v — 3 massive neutrinos, each carrying mass m v . 



Flavors 


Mass limit [eV] 


3.25 


2.95 


3.50 


1.67 


3.75 


1.21 


4.00 


1.02 


4.50 


0.87 


5.00 


0.82 


5.50 


0.77 


6.00 


0.71 



Note that for less than one massive species the mass limit becomes quite poor 
because the limit on Vt v h 2 remains almost constant down to N v = 3.25, whereas the 
limit on m u is proportional to 1/ (N u — 3). 

3.4- Understanding the Q u -N u - degeneracy 

In all of our cases there is a conspicuous degeneracy between Q u and N v . This can 
be understood by studying the power spectra for different models. Figure 4 shows 
power spectra for various ACDM models with = 0.05, Q = 1, h = 0.7, n s = 1, 
and iV^massive = 1, but with different values of Q m , Q u , N u , and b. The full line 
corresponds to the standard ACDM model with no neutrino mass. The dashed line 
is for VL V = 0.05, Q m = 0.25, and = 3, and the power suppression below the free- 
streaming scale is evident. The dotted line is for Q u = 0.05, Q m = 0.25, and N u = 8. 
Here, the spectrum is suppressed even further, at scales up to the horizon size at matter 
radiation equality (because of later matter-radiation equality). Finally the long-dashed 
line is for Q u = 0.05, Q m = 0.35, and N u = 8. Increasing the matter density leads 
to earlier matter domination, which in turn leads to less small-scale suppression and 
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smaller horizon size at matter radiation equality. The power spectrum of this model is 
almost indistinguishable from the standard ACDM model (f2„ = 0, N v = 3), explaining 
the f^-iV^-degeneracy. 



1000 



^ 100 

s; 

10 

1 

0.01 0.10 
k [h/Mpc] 

Figure 4. Power spectra for ACDM models with fib — 0.05, fl = 1, h = 0.7, n s = 1, 
and iV^massivc = 1 and a common large-scale normalization. The full line is for fl u = 0, 
fl m = 0.25, N v = 3, dashed is for fl v = 0.05, fl m = 0.25, N v = 3, dotted is for 
Q v = 0.05, n m = 0.25, N v = 8, and long-dashed is for fl v = 0.05, tt m = 0.35, N„ = 8. 




4. General Thermal Relics 

4-1. Low-Mass Fermions 

Our second generic class of models is that of some new particle with mass m x that 
was once in thermal equilibrium, but decoupled when the universe had a temperature 
T D and the radiation density was characterized by effective thermal degrees of 
freedom. The relevant masses are of order eV so that the particles decouple when 
they are relativistic, i.e. at decoupling they are characterized by a Fermi-Dirac or Bose- 
Einstein distribution of temperature Td. The N v = 3 ordinary neutrinos are now always 
assumed to be massless. 

In the present-day universe, the new particles will be non-relativistic and contribute 
a matter fraction 

2 = ™ | 1 fermions 

183 eV (?* D | | bosons 1 ; 

where gx is the number of the particle's internal degrees of freedom while g w is the 
effective number of thermal degrees of freedom when ordinary neutrinos freeze out with 
g* v = 10.75 in the absence of new particles. For the case of N u ordinary neutrinos with 
equal masses (degenerate mass scenario) we have gx = 2N U (two spin degrees for each 
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flavors), T D »i 2 MeV and thus (^d = 9*u s° that Eq. (4) corresponds to the usual result 
Q u h 2 = N u m u /91.5 eV as stated already in Eq. (1). 

In the late epochs that are important for structure formation, the momentum 
distribution of the new particles is characterized by a thermal distribution with 
temperature Tx that is given by 

T u \9*dJ 

We have modified CMBFAST to incorporate a fermionic thermal relic with such a 
momentum distribution. 

Our first case is gx = 2, i.e. a Majorana fermion with interactions that can be 
weaker than those of ordinary neutrinos. The limits are calculated with a Fermi-Dirac 
distribution, but they would not change significantly for the Bose-Einstein case. In 
Fig. 5 we show the 68% and 95% likelihood contours for the allowed parameters in 
terms of Vt x h 2 and T x /T v . Limits for selected values of T x /T u can be found in Table 5. 
Of course, in the limit Tx/T u = 1 we recover the 3+1 scenario of 3 ordinary massless 
neutrinos and one extra flavor with mass m v . 

Figure 5 reveals that the bound on Qxh 2 is almost independent of Tx down to 
a critical value. For lower Tx the bound abruptly relaxes, allowing for a large mass 
fraction fix of the new particles. This behavior is explained by our non-linearity cut 
in the data at k — 0.2 h Mpc -1 . For low values of T x the free-streaming scale of the 
particles is simply too small to affect large-scale structures above the cut. 

In order to substantiate this interpretation we note that the free-streaming scale 
can be approximated as 

2tt 20Mpc (T x \ l i 

Afs = I7s ~ TirfT I X ) | 1 + log 

In the right panel of Fig. 5 we show contours for which A; FS = 0.1, 0.2, and 0.3 /iMpc~ 
Evidently /cfs = 0.2/iMpc^ 1 corresponds fairly well to the value of T x where the 
constraint on Vtxh 2 disappears. An exact match can not be expected because our 
estimate was only approximate. Moreover, the X-particles follow a thermal distribution 
so that some fraction will have large energies and therefore stream further than the 
mean of the distribution. Therefore, there is some free-streaming effect above our cut 
even if the mean free streaming scale is too small. 

We note that even though the galaxy data from surveys such as 2dFGRS or 
SDSS must be cut at 0.2 h Mpc" 1 , data from the Lyman-a forest measure an earlier 
cosmological epoch and thus the linear regime extends to much larger values of k. 
For instance, this was used in Ref. [29] to constrain the possible mass of a thermally 
produced warm dark matter particle to be mx ^ 750 eV, assuming that it provides all 
the dark matter. In our paper we have deliberately avoided the Lyman-a data since 
the conversion of the measured flux power spectrum into a matter power spectrum is 
fraught with difficulties and the result is at present highly controversial. 
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Figure 5. Likelihood contours (68% and 95%) for a generic thermal relic with gx = 2, 
characterized by its density, flxh 2 , and effective temperature, Tx- The dashed line 
shows A\ 2 = 4, corresponding to the 95% confidence limit for fixed Tx- In the left 
panel, the solid lines correspond to the indicated values of the particle mass in eV, 
whereas in the right panel they denote the parameters for which the free streaming 
scale according to Eq. (6) is kp$ = 0.1, 0.2, and 0.3 h Mpc -1 . 



Table 5. Limits (95% CL) for gx — 2, i.e. a general Majorana fermion. 

Upper limit 



log(T x /T„) 


T x /T u 




n x h 2 


mx [eV] 


0.00 


1.000 


1.00 


0.011 


1.0 


-0.05 


0.891 


1.41 


0.012 


1.6 


-0.10 


0.794 


2.00 


0.015 


3.0 


-0.15 


0.708 


2.82 


0.017 


4.1 


-0.20 


0.631 


3.98 


0.021 


7.5 


-0.25 


0.562 


5.62 


0.029 


15.0 


-0.30 


0.501 


7.94 


No bounds. 



For a particle freezing out after the QCD phase transition, the effective number 
of thermal degrees of freedom at decoupling is g*D^20 or g^j g* u <^2 so that the 
mass bounds are quite restrictive, m x < 3 eV or better. However, at the QCD phase 
transition, g* decreases dramatically due to the large number of confined colored degrees 
of freedom. At epochs shortly before the transition g* > 60 so that a particle freezing out 
at that time is characterized by g^j g* u ^6, preventing us from stating any significant 
bounds without using the Lyman-a data. 

4-2. Scalar Bosons 

The cosmic mass fraction of our new particles is proportional to gx^x so that the limits 
of the previous section can be roughly scaled to cases with a different number of internal 
degrees of freedom gx- Also, on the basis of the estimated free- streaming scale one can 
estimate for which value of Tx/T u the mass limits begin to fail. Still, we have performed 
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Figure 6. Same as Fig. 5 for scalar bosons (gx = !)■ 



Table 6. Limits (95% CL) for gx = 1, i.e. a general scalar or pseudoscalar boson. 

Upper limit 

]oe(T x /T v ) T x /T v g*n/9* v V x h 2 m x [eV] 



0.00 


1.000 


1.00 


0.009 1.2 


-0.05 


0.891 


1.41 


0.011 2.0 


-0.10 


0.794 


2.00 


0.013 3.3 


-0.15 


0.708 


2.82 


0.020 7.8 


-0.20 


0.631 


3.98 


No bounds. 



an explicit analysis for another special case, the one of a scalar or pseudoscalar boson 
{gx = 1)- This example is of particular interest for axions or axion-like particles. 

The results of our likelihood analyis are shown in Fig. 6, which is fully analogous 
to Fig. 5. We also tabulate explicit m x limits for selected values of T x /T u in Table 6. 
We note that for particles decoupling after the QCD phase transition the mass limit for 
equivalent decoupling epochs is only slightly weaker than it is for the previous gx = 2 
case, i.e. the simple mx oc g x X scaling would have yielded less restrictive limits. The free- 
streaming scale enters the non-linear regime at higher Tx than before. Once more we find 
that there is no meaningful mass limit for particles decoupling before the QCD epoch. 

As a concrete example we mention axions. The energy-loss limit from SN 1987A 
implies that the axion mass is bounded by m a <^ 10~ 2 eV [30]. Therefore, it is usually 
assumed that cosmologically significant axions must be non-thermal relics with very 
small masses. On the other hand, the SN 1987A argument may have loop holes or one 
may assume that the axion couplings, notably those to photons, are quite different from 
the most common models, perhaps allowing eV-mass axions to exist in the "hadronic 
axion window" [31]. If the axion mass is in the eV range, its couplings to nucleons and 
pions is so strong that it decouples thermally well after the QCD transition [32-34]. 
According to Table 6 this implies a cosmological mass limit of m a < 2-3 eV. 

In a direct telescope search, narrow lines from the two-photon decay of axions 
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trapped in galaxy clusters were searched in the mass range 3-8 eV [35], but no signal 
was found. According to our new limits, axions in this mass range can not exist because 
of their excessive free-streaming effect on cosmic structures. 

On the other hand, the possibility of axion-like particles that couple only to photons 
by virtue of their 7r°-like two-photon vertex is sometimes discussed [36,37]. For an 
interaction strength small enough to obey the well-known globular-cluster limit [30], 
such particles would decouple well before the QCD epoch [36]. Therefore, our arguments 
do not provide any new limits on such hypothetical particles. 

Either way, our limits do not infringe on the possibility that the CAST experiment 
at CERN to search for solar axions could discover axions or axion-like particles with 
masses up to the eV range [38]. Actually, a small gap will remain between our 
new limits and CAST's future high-mass frontier. Likewise, the future possibility to 
search for axion-like particles with a string of decommissioned HERA magnets [39] 
remains unaffected. 

5. Discussion and Summary 

We have calculated general cosmological bounds on the masses and abundances of light, 
stable particles such as neutrinos and axions. For neutrinos we have calculated three 
different generic cases: 

1) All neutrinos are massive, so that there are N v neutrinos, each with mass m v . 

2) Only one neutrino is massive, but altogether there are N v neutrinos. 

3) There are three massless neutrinos and N v — 3 species with mass m v . 

In each case it is found that there is a strong correlation between Vt u h 2 and N v , but not 
necessarily between m v and N v . 

We have extended the standard analysis to include generic light particles such as 
axions and majorons. We assume such particles to have been thermally produced so 
that they are characterized by an effective temperature Tx- The other free parameter 
can be taken to be either VLxh 2 or mx, but since Vtxh 2 is the main parameter entering 
into CMB and LSS calculations we chose the parameters to be flxh 2 and Tx- 

In our analysis we do not use LSS data in the non-linear regime, i.e. at small scales, 
and we also do not use Lyman-a data. Therefore, particles that decouple too early 
and thus today have much smaller momenta than neutrinos do not produce significant 
free-streaming effects on the relevant scales. For both spin-| fermions and for scalars 
that decouple before the QCD epoch, we can not state a significant mass limit. On the 
other hand, for such particles decoupling after the QCD epoch, there is a mass limit 
between 1-3 eV, depending on the spin and exact epoch of decoupling. For axions, our 
limits supersede previous results from direct telescope searches for axion decays. 

An independent study of cosmological bounds on neutrino masses and relativistic 
relics was completed shortly after ours [40]. For overlapping cases there is agreement 
with our work. 
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